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ABSTRACT 



An Advanced Simulation of Digital Filters has been 
implemented on the IBM 360/67 computer utilizing Tektronix 
hardware and software. The program package is appropriate 
for use by persons beginning their study of digital signal 
processing or for filter analysis. The ASDF programs provide 
the user with an interactive method by which filter pole and 
zero locations can be manipulated. Grpahical output on both 
the Tektronix graphics screen and the Versatec plotter are 
provided to observe the effects of pole-zero movement. 
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I. INTRODUCTION 



Digital methods have recently begun to dominate the 
selection of transmission and filtering schemes. When the 
student first confronts the discrete domain, he encounters a 
series of unfamiliar concepts. While most of the discrete 
concepts have parallels in the continuous or s-domain, this 
dualism Is not apparent initially. 

The programs comprising the Advanced Simulation of 
Digital Filters (ASDF) provide the user with a convenient 
method for developing intuition In the discrete domain. 
This is achieved by allowing the user to generate digital 
filters by directly entering the ooles and zeros into a 
graphics terminal. The ASDF programs then perform the 
calculations required to display the cha racte r i s t i cs of the 
user generated filter. 

The overriding objective of the ASDF program package is 
to provide the user with an interactive method for modifying 
the filter pole/zero locations in the z-plane, as well as to 
compute the filter responses. Both time and frequency 
responses are computed and displayed graphically on the 
Tektronix 4012 graphics terminal. [Ref. 9] 
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II. SELECTION OF ASDF ALGORITHMS 



There are several parameters which are particularly 
important in selecting the algorithms used in the ASDF 
programs. Since one of the primary objectives of these 
programs is to provide the user with clear explanations of 
methods for computing digital filter characteristics, direct 
or straightforward methods should be used. The importance 
of developing a user's intuition for the discrete domain and 
digital filters suggests that graphical rather than 
numerical outputs would be most appropriate. The use of 
graphical outputs to characterize a filter's responses 
implies minimum importance of precision in these results. 
Since the reason for computing the filter characteristics is 
to develop the user's Intuition for the discrete domain, 
there are no pressing constraints on the selected algorithms 
with regard to computational efficiency, speed, or memory 
space utilization. In fact, the final ASDF package spends 
the greatest portion of its execution time making the 
interface with the z-plane, and the form of the output easy 
to understand and to use. The greatest part of the memory 
space utilized for the ASDF programs is consumed in the 
interfacing processes rather than in the methods by which 
the data is stored, used, or computed. 
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A. ALGORITHMS SELECTED 

Because of the above stated objectives, a simple 
algorithm was selected for computation of time domain 
responses of digital filters. Once the user enters the 
desired pole and zero locations, a dl rect-form-two 
difference equation representation of the filter is 
generated [Ref. 5]. The unit sample sequence is used for 
the filter input, the difference equations are solved 
recursively In the discrete time domain and the unit sample 
response is recorded as the output. The unit step response 
is calculated by a similar method. 

It is known that the di rect-form- two filter 
implementation is very sensitive to parameter variations, 
and as such cannot be utilized in situations requiring 
extreme accuracy. The advantages associated with simplicity 
of understanding the algorithms, the ease of programming the 
general case, and the deemphasis of accuracy in outputs for 
marginally stable filters outweighs this disadvantage. 

B. DEFINING THE DATA STRUCTURE 

A detailed description of the important variables in the 
ASDF programs is given in Appendix B. There is, however, an 
overlding data structure common to all portions of the ASDF 
package which is important. 

The array POLZRO holds the locations of the poles and 
zeros of the filter under analysis. In all cases, zeros are 
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stored in P0LZR0(l,J) for values of I from one to ten. Pole 
locations are stored in POLZRO(I,J) for values of I from 
eleven to twenty. Only the top half plane pole of each 
complex pair is stored. There are five sections of 
POLZRO(l/J)/ i.e. values of J from one to five. POLZRO(I^l) 
is the real part of- the rectangular representation of the 
i th root. P0LZR0(I,2) is the imaginary part of the 
rectangular representat ion of the i th root. P0LZR0(I,3) and 
P0LZR0(I,4) are the magnitude and angle portions of the i th 
root location in polar form. P0LZR0(I^5) is the root type 
designator. Root type designators are: one for real zeroS/ 
two for complex zeros, three for real poles, and four for 
complex poles. 

The array IROOTS(I) holds the numbers of roots of each 
type of root already entered into the system. IROOTS(l) and 
IR00TS(3) hold the number of real zeros and poles, 
respectively. IR00TS(2) and IR00TS(4) hold the number of 
complex zero and pole pairs, respectively. 

The variable I ERROR is the error flag, which is set 
whenever the user generates an error. There are numerous 
possible errors. These include trying to delete roots which 
have not been entered, or adding roots which cause the 
maximum system order to be exceeded. 
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IMPLEMENTED SOLUTION 



Having determined the general form of the solution it 
remained to generate the software implementation. There 
V'/ere five major areas to be programmed: the interface with 
the z-plane, the time response, the frequency response, the 
plotting of output curves, and the supervisory routines. 

A. INTERFACE WITH THE Z-PLANE 

Since the method by which the user modifies the pole and 
zero locations in the z-plane is of no theoretical interest, 
the discussion of this implementation will be brief. 
Extensive coverage of the subroutines employed and the 
definitions of the variables used is reserved for Appendix 

B. 

The ASDF command which allows the user to make pole-zero 
modifications in the z-plane is POLZRO. This command 
initiates execution of a load module which is structured as 
depicted in figure 3-1. Each of the boxes shown reoresents a 
subroutine. Lines with arrow heads on both ends indicate 
that program control passes to the subroutine and then 
directly back to the calling routine with no intermediate 
branches. The main program, MAIN, initializes the required 
arrays, and tests to ensure that the user has generated a 
realizable causal filter. 
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START 




Figure 3-1 Flowchart Of POLZRO 
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The actual Interaction with the z-plane is accomplished 
in subroutine UMTCIR. The structure of UNTCIR is diagrammed 
in figure 3-2. The large double ended arrow indicates that 
there are many subroutines of the form CRRMMM. The 

subroutines LETTER and BOXUC generate the unit circle and 
command array. PLOT-10 software [Ref. 10] statements 

determine the location of the user positioned cursor. This 
cursor position is adjusted with software tabs to one of 
twenty-six possible locations. Based on the X and Y 
location of the adjusted cursor position a specific 

subroutine is selected for execution. The flashing 
alphanumeric cursor is positioned in the selected command 
box using the adjusted cursor location as a reference. The 
subroutine then reads two alphanumeric characters typed by 
the user. The combination "space”, "e" indicates that the 
user wishes the selected subroutine to be executed. The 
combination "space", "x" indicates that the user wishes to 
terminate the pole zero manipulation phase of the program. 

All of the commands which change the status of the 
z-plane have names of the form CRRMMM where C is either "A" 
for add root(s) or "D" for delete root(s). The RR indicates 
root type; RZ - real zero, CZ - complex zero, RP - real 
pole, CP - complex pole. The MMM stands for the mode of 
data entry. lAG stands for interactive graphics i.e. the 
roots v/i 1 1 be adjusted with the cursor. POL and RCT stand 
for polar and rectangular, respectively. The general 
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Figure 3-2 Structure Of Subroutine UNTCIR 
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structure of the CRRMMM subroutines is described in figure 

3 - 3 . 

All CRRMMM subroutines follow the same general pattern. 
The POLZRO table is checked to determine whether or not the 
issued command is valid. Poles and zeros must be present 
before they can be deleted. There cannot be more than ten 
poles or ten zeros. The desired pole or zero location is 
next acquired. For subroutines adding roots^ the entered 
location is that of the root(s) to be added. For 

subroutines deleting roots the entered location is any 

location near the root(s) to be deleted. Using subroutine 

LOCATE, the POLZRO table is searched for the root of the 

correct type which is nearest the snecified location. This 
nearest root is then highlighted. In all cases, except when 
adding roots in the I AG mode, the user is asked whether or 
not the correct location has been determined. Once the root 
location is verified the root is added and plotted with 
subroutines STRRTS and PLOTRT. For subroutines deleting 
roots the POLZRO table is adjusted with subroutine UPDATE 
and the unit circle is redrawn. 

The user may make an error in carrying out the above 
steps, i.e., by placing a pole outside the unit circle, by 
exceeding the maximum system order, or by entering a zero 
outside ten units from the origin. In each of these cases an 

^Hereafter, poles and zeros will be refered to as roots. 
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Figure 3- 3a Structure of CRRMMM Subroutines 
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Figure 3-3b Structure Of CRRMMM Subroutines Continued 
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Figure 3- 3c Structure Of CRRMMM Subroutines Continued 
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error fla» is set. Subroutine ERROR handles the printing of 
messages and subroutine 1ST prints a listing of the current 
filter status before the program resumes. 

B. TIME RESPONSE 

The transfer function of a digital filter is defined in 

equation 3.1. The z and b^ will be defined in a 

subsequent section. The variable A is an arbitrary constant 

assigned by the user and is called the filter gain. 

M M ,, 

A IT ( z-z . ) Zb z^"^ 

i=l ^ r=0 

H(z) => = (3.1) 

N N _ 

IT (z-p.) 1- Z 3 Z~ 

i=l ^ k=l ^ 

Inherent in generating the time domain response with a 
di rect-form-two filter implementation are the determi nat ion 

of the coefficients a, and b in equation (3-1), and the 

k r 

computation of the unit sample and unit step sequences. 
Having determined the filter coeficients, the time domain 
response of the filter is determined by using the 
appropriate filter input sequence and recording the output 
of the filter for the desired number of points. This is 
accomplished by recursively solving a difference equation 
which character i zes the filter for 1024 time points. The 
standard di rect-form- two filter implementation is shown in 
figure 3-4. The z~^ terms represent unit time delays. This 
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Figure 3-4 Direct Form Two Filter Implementation 






figure corresponds to the difference equations: 



N N 

y(n) = I a, y(n-k) + Z b x(n-N+r) (3.2) 

k=l ^ r=0 

To compute the filter's time response, three operations 
are required: calculate the and coefficients, generate 
the input sequence x(k), k=0, 1, . . . 1024, solves the 

difference equations and stores the filter output values. 
There are four central subroutines which perform these 
operations: ASUBK, BSUBR, RSPNSE, and COEFF. Subroutines 

ABSUK, BSUBR, and COEFF generate the filter coefficients 
starting with the pole and zero locations specified. 
Subroutine RSPNSE solves the difference equations by 
coordinating the other three routines, generating the 
appropriate input sequence, and storing the output sequence 
val ues . 

The ASUBK and BSUBR subroutines are identical in 
structure and purpose. ASUBK computes the degree of the 
denominator polynomial as two times the number of complex 
pole pairs plus the number of real poles. A complex vector 
is formed named POLES, which holds a listing of all the 
poles of the filter under analysis. The POLZRO array stores 
real poles in the form X+jO. Only the top half plane pole 
of each complex pole pair is stored, i.e., (X+jY), where X 
and Y are non-negative. The vector POLES has entries for 
all real poles and both poles of each complex pole pair. 
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Once the POLES vector has been generated the subroutine 
COEFF multiplies out the factors of the form (z-p^) to 
generate the coefficients of the denominator polynomial. 

The algorithm for multiplying out the denominator 
polynomial is essentially the same as would be done by hand, 
multiplication of a series of literal factors of the form 
(z-P^). The vector POLES is passed to the COEFF subroutine 
as the vector V. Since the denominator factors are of the 
form (z-p^)/ the V vector is changed to hold the values • 
If the system is of order less than ten, the remaining 
positions of V are filled with zeros. The degree of the 
denominator polynomial Is subsequently used to determine how 
many poles are actually located at the origin. 



The vectors 


used 


1 n 


the subrout i ne 


COEFF 


hold the 


coefficients for 


the 


po 1 ynom i al 


bei ng 


generated 


The 


literal powers of 




are 


supressed 


. The 


vector 


HI 


always 



contains the partial product of the K factors previously 
multiplied together. The vector H2 represents the old 
partial product times the "z” of the (z-Pi> factor. Since 
the coefficient of "z" is always one, H2 is a left shifted 
version of HI. The vector HS represents the partial product 
times the next pole location (Pj^). This next -p^ value is 
stored in V(l). The new HI vector is formed by adding the 
H2 and H3 vectors. HI now contains the product of the first 
K+1 factors. The values in the V vector are next shifted 
such that the next pole location to be multiplied is always 
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mm 








in the V(l) position. These four steps are repeated ten 
times to j^enerate the coefficients of the complete 
polynomial in "z". The formula for the transfer function 
requires coefficients of the denominator polynomial in terms 
of z~^ so the order of the coefficients is reversed. 

Figure 3-5 depicts several iterations of the algorithm 
used in the COEFF subroutine. Figure 3-5a shows the initial 
conditions for each of the vectors defined in the above 
paragraph. Figure 3-5c shows the values of the vectors 
after the first iteration. Figure 3-5f shows the vector 
values after the second iteration. The third and subsequent 
iterations follow this pattern. 

In precisely the same way as di scribed above, the a^^ 
coefficients are generated with the ASUBK subroutine. The 
zero locations are stored in POLZRO exactly as are the 
poles. The generation of the vector ZEROS and the use of 
the subroutine COEFF parallels the BSUBR subroutine exactly. 

The significant variables in subroutine RSPNSE are: 
XSUBN, INPUT, OUTPUT, FLTRGN, A, and B. The vectors A and B 
hold the already computed values of aj^ and bj^. A(l) holds 
ag; A(2) holds a^^ . . . A(ll) holds The b^^ are stored 
similarly in the B vector. XSUBN is the current value of the 
input sequence. INPUT is the value of the node so labeled 
in figure 3-4. OUTPUT is similarly noted on figure 3-4. 
FLTRGN is the value of the filter gain, and is assigned by 
the user. The vector Z holds the delayed INPUT values. 
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Figure 3-5a. Initial Conditions Of COEFF Vectors. 




Figure 3-5b. Values Of COEFF Vectors Before The End Of The 



First Iteration. 
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Figure 3-5c. Values Of COEFF Vectors starting The Second 
I teration. 
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Figure 3-5d. Values Of COEFF Vectors before The End Of The 
Second I teration. 
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Figure 3 - 5 e. Values Of COEFF Vectors Starting The 
Third I teration. 
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Figure 3 - 5 f. Values Of COEFF Vectors before The End Of The 
Third I tera t i on . 



26 



Z(l) is no delay; Z(2) is delayed one time unit; Z(3) is 
delayed two time unitS/ etc. 

The subroutine RSPNSE generates the appropriate input 
sequence. For the unit sample response the input sequence 
is 1, 0, 0/ 0^ ... . For the unit step response the input 
sequence is 1, 1, 1, ... . The appropriate value of 

XSUBN is determined by the type of input sequence. The 
variable INPUT is then computed to be the sum of XSUBN and 
the a.'s times the delayed INPUT values. The filter output 

iC 

is the sum of the b 's times the delayed INPUT values. 

r 

These output values are then multiplied by the filter gain, 
FLTRGN, and stored in the appropriate array: Yl for the 
unit sample response, YS for the unit step response. 

C. FREQUENCY RESPONSE 

For discrete systems, the frequency response is found by 
evaluating the system transfer function at points on the 
unit circle. Subroutine FRQNCY generates both the phase 
and the magnitude of the transfer function. 

The important variables in the FRQNCY subroutine are: 
AA, BB, THETA, Z, NNUM, and NDEN. The variables AA and BB 
are the double precision values of the A and B vectors 
described in the section on time response. THETA is the 
angular displacement around the unit circle for which a 
value of the frequency response is calculated. Z is the 
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classical z-plane variable: 



Z=exp(j THETA) (3.3) 

NNUM is the number of numerator terms. NDEN is the number 
of denominator terms. 

The frequency response Is computed for 1024 points 
around the unit circle. All points computed are for THETA 
{greater than zero and less than tt. The values of Z for 
THETA greater than -tt and less than zero are neither 
computed nor plotted since the magnitude of the transfer 
function is an even function of THETA and the phase response 
is an odd function of THETA. 

The generation of the a^^ and b^ coefficients has already 
been discussed in the section on time domain response. The 
evaluation of H(z) is straightforward. The value of THETA is 
incremented through 1024 steps from zero to tt . For each 
value of THETA generated, the value of Z is computed. The 

value of H(z) is then computed in two parts. The numerator 

— p 

is the sum of NNUM terms of the form b^ times z for 

r=0, 1, 2,. . . NNUM-1 \ The denominator is one minus the 

sum of (NDEN-1) terms of the form a^^ times z for 

2 

k=l, 2, 3, . . . NDEN-1 . Dividing numerator by denominator 

^ NNUM Is equal to the degree of the numerator plus one. 

2 

NDEN equals the degree of the denominator plus one. 
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yields H(z). The phase of H(z) is then determined using the 
principal values of the arctangent function; the results are 
stored in the vector YP. The magnitude of H(z) is computed, 
multiplied by FLTRGN, the filter gain, and subsequently 
stored in the vector YM. 

D. OUTPUT PLOTS 

The output plots generated by ASDF are generated by two 
different software programs each being different 
implementations of the same plotting requirements. One set 
of plots is generated using PLOT-10 software on the 
Tektronix 4012 terminal. Another set of plots utilizes the 
Versatec plotter and associated software. [Ref. 7] The 
PLOT-10 software is adequately covered in reference 10. The 
Versatec software is explained in Technical Note-034 in the 
computer center. 

The plots which are generated on both devices are 
essentially the same. There are five plots generated on the 
4012 terminal: unit sample response, unit step response, 
phase response, magnitude of the transfer function, and 
magnitude of the transfer function in decibels. Both of the 
time domain plots are terminated when the output value has 
stabilized to within one percent of the maximum generated 
system output. All 1024 points of the phase and magnitude 
plots are plotted. 
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The form of the Versatec output is similar to that of 
the 4012. One additional plot is included. It shows the unit 
circle with poles and zeros of the filter being analyzed. 
Plotting is terminated when the filter output stabilizes to 
within one percent of the maximum system output generated. 

On both devices, the values of unit sample response are 
stored in a vector YI and plotted with subroutine PLTIMP. 
The unit step response values are stored in YS, and plotted 
with subroutine PLTSTP. The phase response is stored in 
vector YP. The magnitude of the transfer function is stored 
in vector YM. The magnitude of the transfer function in 
decibels is stored in YMDB. The three vectors YP, YM, and 
YMDB are all plotted with subroutine PLTTRF. The Versatec 
plot of the unit circle is generated with subroutine UNTCR. 
The 4012 plot of the unit circle and command array are 
generated with the subroutine UMTCIR. 

E. SUPERVISORY ROUTINES 

There are three supervisory routines used with the ASDF 
programs. These supervisory routines make the loading of 
ASDF transparent to the user. There is a substantial amount 
of file manipulation required to store and load the ASDF 
programs. A knowledge of the hardware and software 
available for these purposes is not oar t i cul ar 1 y useful to 
the user. For these reasons, the supervisory routines are 
discussed only briefly. Two of these routines are based on 
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CP/CfIS and IBM EXEC environment; the third is written in the 
IBM-360 Job Control Language. 

The JCL program moves the ASDF program from its 
permanent storage location on DISK02 to the CP/CMS disks 
OSXOOl. The two EXEC routines generate an appropriate 
CP/CMS work space, and manipulate the ASDF text file into a 
usable form. 

The EXEC program ASDF191 resides on the user 191-P-disk. 
This routine automatically requests and formats ten 
cylinders of temporary disk space. The ASDF TEXT file is 
moved from the OSXOOl disk to the newly created temporary 
disk space 192-T. Once in place. The ASDF TEXT file name is 
changed to ASDF EXEC to generate the next file with the 
correct file type. ASDF EXEC is then edited to remove the 
first set of data cards which hold the ASDF192 EXEC program. 
The 192-T-disk is released and the user is requested to type 
two commands to the system. The "login 192 P" command makes 
the 192 disk the P-disk. The "EXEC ASDF192" command begins 
the execution of the second EXEC file. 

The general ASDF192 EXEC structure is depicted in figure 
3-6. ASDF192 EXEC is created by ASDF191 EXEC and resides on 
disk 192. In a manner similar to ASDF191, ASD'^192 splits 
the ASDF TEXT file, which contains all of the component 
parts of the ASDF program, into a series of smaller files. 
Once all files have been separated, titles are altered 
appropr i atel y, text files are loaded into core, and load 
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Figure 3-6 Structure Of ASDF192 EXEC 
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modules are generated. ASDF192 then reverts to an infinite 
loop/ reading and implementing the desired user command. 
There are six user commands available to the user: POIZRO, 
RESPONSE/ EXEC PUNCH/ EXEC HRD$CPY/ OPTIONS/ and INTRO. The 
applicability of these commands is briefly discussed in 
chapter four. This brief description is also presented by 
executing the OPTIONS command once ASDF has been loaded. 
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IV. STRUCTURE OF ASDF COMMANDS 



The ASDF program generates the time and frequency 
responses of realizable digital filters of order ten or 
less. Chapter three describes the algorithms used for these 
computations. The procedure for executing the ASDF program 
generally consists of tv/o distinct steps. The first step is 
to enter the poles and zeros of the filter to be analyzed. 
The second step is to perform the calculations required to 
generate and plot the desired filter responses. There are 
six commands which may be executed v/ithin the ASDF 
environment providing the user with the ability to control 
the form of the output, to enter the poles and zeros of the 
filter, and to clarify most questions pertaining to the 
meaning and use of the general commands. The six 
fundamental commands available to the user are: OPTIONS, 
INTRO, POLZRO, RESPONSE, EXEC PUNCH, and EXEC HRDSCPY. Only 
the POLZRO command has a subordinate command structure. All 
of these commands are issued by typing the command on the 
terminal keyboard in response to: *ASDF-RE. 

The OPTIONS command provides a listing of all commands 
available to .the user. Each command is accompanied by a 
brief description of its intended use. 

The INTRO command v/ri tes on the screen .six pages of 
information pertaining to the POLZRO command substructure. 
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INTRO includes a general description of the data input 
nodes, and a description of several input situations with 
the appropriate executable commands. 

The POLZRO command allows the user to modify or create a 
filter in the z-plane. All adjustment of poles and zeros 
must take place with this command. By Issuing this command 
the user may add and delete poles and/or zeros as well as 
adjust the filter gain constant. The commands which are 
issued within the POLZRO environment are described in a 
later section of this chapter. 

The RESPONSE command displays the character i s t i cs of the 
filter being created by the user. Five displays are 
presented: unit sample response, unit step response, phase 
of the transfer function, magnitude of the transfer 
function, and the magnitude of the transfer function in 
dec 1 be 1 s . 

The EXEC PUNCH command allows the user to punch a card 
deck of the filter pole and zero locations on the offline 
punch. This punched output can be nicked up in room 1-140 
of the computer center. The cards punched by this command 
may be subsequently reentered into the system so that this 
particular filter can be referenced or modified at a later 
date. 

When the user desires to enter a filter in this way, the 
punched cards are presented to one of the computer operators 
over the counter in room 1-140. The cards are punched in 
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the correct format for direct submission into the CP/CMS 
card reader. The operator subsequently reads this deck into 
the user's virtual card reader. As the user signs onto a 
terminal and gains access to CP/CMS/ the computer includes 
an additional message in the standard login typeout showing 
that the user has a virtual reader file. This message is of 
the form: 

FILES:- 01 RDR, NO PRT, NO PUN 

When the user wishes to access this old filter data he must 
execute the CP/CMS command OFFLINE READ * The OFFLINE 
READ command will read the punched card filter data from the 
user's virtual card reader, and generate a file on the 
user's disk space entitled FILE FTOIFOOI. Any files with 
this name present before the OFFLINE READ command is issued 
will be ovei — written by the new file. If the user does not 
v/ish to lose an older FILE FTOIFOOI, he should use the 
CP/CMS ALTER command to change the name of the old file that 
he wishes to protect. Once the user has successfully loaded 
the ASDF program, the ASDF signature i.e., *ASDF-RE, is 
presented. The OFFLINE READ command may be issued any time 
that this signature appears. 

The EXEC HRD;?CPY command is available only when the 
SUBMIT function of CP/CMS is functioning. When the SUBMIT 
command is not functional a message is sent to the user 
during the login procedure v/hich states that SUBMIT is not 
available. the output from this command is ultimately 
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available in the computer center, and consists of two parts. 
The first part of the output is the numerical listing of the 
values of the responses computed by the ASDF program. The 
EXEC HRD$CPY command is the only command which will allow 
the user access to these numerical values. The second part 
of the output is a series of Versatec plots, each similar to 
the screen presentations of the RESPONSE command, with one 
additional plot of the unit circle and pole/zero locations. 
When the EXEC HRD5CPY command is executed the ASDF program 
will ask several questions, using the responses to generate 
a JOB card. The response to the ’’ENTER DESIRED JOB NAME” is 
the title which appears on the final output. ASDF then 
submits the filter character i st i cs to the operating system 
for processing. The hard copy output v'/ill be found in the 
standard output boxes. 

A. SUMMARY OF DATA INPUT MODES 

There are three nodes through v/hich the poles and zeros 
may be entered into the z-plane: 1) interactive graphics, 
2) rectangular, and 3) polar. Each of these modes has 
particular attributes and limitations. 

1. Interactive Graphics (lAG) 

This mode was designed to enter poles and zeros 
rapidly, in cases where accuracy Is not of prime importance. 
The interactive graphics mode is characterized by the 
appearance of the graphics cursor (the intersection of two 
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lines which trace across the terminal screen). The position 
of the cursor is controlled by the two thumb wheels located 
to the right of the interactive ASCII keyboard on the 
Tektronix 4012 terminal. The coverage of the interactive 
graphics cursor is limited to values of X and Y greater than 
-1.1 and less than +1.1. Zeros outside this area must be 
entered with either the polar or rectangular commands. 
Careful adjustment of the cursor can usually locate poles or 
zeros within .002 units of the desired location. Poles must 
be located within the unit circle. 

2 . Rectangular 

This mode was designed to enter from the terminal 
keyboard, poles and zeros which are known in rectangular 
coordinates. It is slower than the interactive graphics 
mode because there is a significant amount of screen 
rewriting required. The rectangular mode allows much 
greater accuracy in the location of ooles and zeros. Poles 
may not be entered on or outside the unit circle. Zeros must 
be located within ten units of the z-plane origin. 

3. Polar 

This mode was designed to enter from the terminal 
keyboard, poles and zeros which are known in polar 
coordinates. This mode is also slower than the interactive 
graphics mode. It allows poles and zeros to be entered to 
any permissible position in the z-plane. Poles must be 
entered within the unit circle. Zeros must be entered within 
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ten units of the z-plane origin. Each polar command has a 
description of the constraints on the data to be entered 
automatically written on the screen. This mode is 
part i cular i 1 y appropriate when constructing all-pass linear 
phase filters. 

B. THE POLZRO COMMAND SUBSTRUCTURE 

In an attempt to maximize execution sneed, a minimum of 
information is "/ri tten on the screen, thus reducing the 
rewrite time after execution of each command. When program 
execution begins, the user is presented with the unit 
circle, and a listing of commands on the top left side of 
the screen. The commands are abbreviated in the following 
manner. The general command type or mode is specified 
across the top of the array as INTERACTIVE, RECT, and POLAR 
(see figure 4-1). Each of these modes conforms to the above 
description. Down the left side of the array are the 
commands of the form CRR, where C is the change to be made, 
A - for add a root, D - for delete a root. RR is for the 
type of root to be acted upon by the command: RZ - real 
zero, CZ - complex zero, RP - real pole, and CP - complex 
pole. As an example, the box in the center column aligned 
with the row titled DCZ stands for the command "Delete a 
Complex Zero - Rectangular Coordinates. There are two 
unique commands listed as MUC, and LST. The New Unit Circle 
(NUC) command simply redraws the unit circle and command 
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UNIT CIRCLE COMMfiNDS 
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Figure 4-1 Basic Display Format 



array. This is appropriate when either too many entries 
have been made to the terminal, and the command array is no 
longer clear, or when noise in the line from the computer 
causes a deformed array or unit circle to be drawn. The LST 
command provides the user with a listing of the current 
status of the filter being analyzed, in a self explanatory 
format . 

In all cases when the unit circle and command array Just 
described appear on the screen, an audible tone or bell is 
heard. This bell indicates to the user that the cursor is 
to be positioned within one of the boxes in the command 
array. In order to execute a command, the cursor is 
positioned in the box associated with the desired command. 
A single keyboard character (any character but the "return" 
key is acceptable) is then typed. The cursor then 
disappears. The system presents a carat oointing to the 
command that the computer is about to execute; the flashing 
cursor blinks within the box associated with the command 
which is about to be executed. If the box so indicated is 
the command that was selected, the user types a "space", an 
"e" (for execute) and then "carriage return". In cases 
where a rectangular or polar mode command is to be executed, 
the screen is then erased, and new information is printed. 
In cases where an interactive graphics command i s to be 
executed, only the cursor reappears. This second appearance 
of the cursor is not accompanied by a bell or audible tone. 
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This indicates that the cursor is to be positioned within 
the unit circle. Once so positioned, a single keyboard 
character is typed (not "return"), and, as appropriate, the 
unit circle is redrawn to reflect the change in pole or zero 
locations. After this update of the z-plane, the circle and 
command array are redrawn with the cursor accompanied by the 
audible tone. 

There are only two exceptions to this scenario. First, 
if the flashing cursor does not flash in the correct box, 
the wrong command will be executed. If an Incorrect box 
contains the flashing cursor, the user should not type the 
"space", "e", sequence. Instead, he should type tv/o other 
characters (e.g., "space", "space") and then "return." The 
cursor returns along with the audible tone, calling for a 
readjustment of the cursor, closer to the desired command 
box. The second exception occurs when the placement of all 
of the poles and zeros is completed. In this case, after 
aligning the cursor with any box, the sequence "space", "x" 
(for exit), "return" should be typed. This will terminate 
the input, or modification, stage of the program and present 
the user with a series of options for the next stage of 
filter processing. At this point a display appears which 
explains the options available to the user. 



C. SUMMARY OF POLZRO SUBCOMMANDS 

Vihen the cursor appears along with an audible tone, the 
user is to position the cursor within the command array, and 
type a keyboard character (not "return”). V/hen the flashing 
cursor appears within a box, he types ("space", "e", 
"return") to execute the specified command. "Space" and any 
other character but "x" will return the user to the command 
table for cursor readjustment. To terminate the input, or 
modification mode, the cursor is positioned over any box, 
and once the flashing cursor enters the box, the characters 
"space", "x", "return" are typed. 
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V. LOADING THE ASDF PACKAGE 



The procedure for using the ASDF program consists of 
either six or seven steps depending on whether or not the 
ASDF package is resident on OSXOOl. Usually the program 
sections of ASDF are stored on the T2314 disk pack OSXOOl. 
Each Monday morning the OSXOOl disks are scratched so that 
the ASDF program must be moved from the permanent storage 
location on DISK02, to the temporary storage location on 
OSXOOl. The user may simply sign onto a 4012 terminal/ and 
execute three commands to load ASDF into his virtual 
machine. Once this loading process has been successfully 
completed, the commands previously discussed are at the 
user's disposal for filter analysis. 

A. LOGIN PROCEDURE 

With one minor exception, the standard login procedures 
(specified in Section 4.4.3 of the CP/CMS User's Manual) 
[Ref. 4] are applicable when entering the CP/CMS system to 
use ASDF. Due to the size of the ASDF programs the user 
must request a 400k virtual machine v/hen loging onto CP/CMS. 
The virtual machine is requested as follows: 

login XXXXTNN 400k 

where XXXX is the user's computer center number, T is the 
type of user i.e. G for general user, and P for private 
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user, and NN Is the two dl^lt terminal number. The 400k 
request in the LOGIN is absolutely mandatory. V/i thout this 
request the user v/ill be able to load the ASDF pro,^ram, 
however, when attempting to execute the RESPONSE command, an 
error: "INPUT COMMAND ERROR" will result. Should this error 
occur when the command is corrctly spelled the user has no 
alternative but to logoff and log back on to the system 
requesting the correct size virtual machine. The remainder 
of the login procedure is identical to the examples in the 
CP/CMS user's manual. 

B. TESTING OSXOOl FOR ASDF PROGRAMS 

The second step In loading the ASDF programs is to 
determine whether or not a copy of the programs exists on 
OSXOOl. The ASDF programs are permanently stored on the 
operating system DISK02. The first time that the programs 
are used each week, the user must move them from DISK02 to 
the disk OSXOOl. To determine whether or not OSXOOl has the 
ASDF package, issue the command: 

TOCP ASDF191 EXEC PI 
The system will respond: 

EXECUTION BEGINS... 

237 = OSXOOl 
ENTER DSNAME OR '/*' 

The correct user response is: 

F0718.asdfl91 
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At this time one of two responses will appear. The 
particular response received is dependent on whether or not 
ASDF resides on OSXOOl. 

1. ASDF is on OSXOOl 

When the ASDF programs are located on OSXOOl, the 
system response is: 

EXTENTS EXHAUSTED, EOF SIMULATED. 

4 BLOCKS 63 RECORDS READ 

This message indicates that file ASDF191 has 
sati sfactor ial 1 y been transferred to the user’s virtual 
machine and the user is ready to continue with the 
subsequent loading steps. 

2. ASDF is not on OSXOOl 

'When the ASDF programs are not on OSXOOl, the system 
response will be : 

DATA SET NOT FOUND (OS/360 S213) 

ENTER DSNAME OR '/*' 

This response indicates that the user must move the ASDF 
programs to OSXOOl before the loading procedure nay be 
resumed. The simplest method for accomplishing this 
transfer is to generate a file on CP/CMS which exactly 
duplicates the Job Control Language program which follows. 

To generate the file which transfers ASDF, logon to 
CP/CMS. Once the system's responses have been entered and 
the system is ready, issue the command: 

EDIT ASDF JCL 
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Then type in, or read 



In , the program which follows. 



// Standard JOB Card (User's Manual Sec. 3. 3. 2.1) 

// EXEC PGM=IE3GENER 
//SYSPRINT DD SYSOUT=A 
//SYSIM DD DUMMY 

//SYSUTl DD DSN = F0718. ASDF,VOL=SER=DI SK02,IJNIT=3330, 

// DCB=(RECFM=FB,LRECL=80, BLKSI ZE=3200 ) , D I SP=(OLD, KEEP) 
//SYSUT2 DD DSN = F0718. ASDF, VOL=SER=OSX001,IJNIT = 2314, 

// DCB=(RECFM = FB, LRECL = 80, BLKS I ZE = 3200 ) , 0 I SP= ( NEW, KEEP) , 

// SPACE=(TRK, (40, 1),RLSE) 

// EXEC PGM=IEBGENER 
//SYSPRINT DD SYS0UT=A 
//SYSIN DD DUMMY 

//SYSUTl DD DSN=F0718. ASDF191, V0L=SER=DI SK02,UNIT=3330, 

// DCB=(RECFM=FB, LRECL=80, BLKS I ZE=800 ) , D I SP=(OLD, KEEP) 
//SYSUT2 DD DSN = F0718.ASDF191, VOL=SER=OSX001,IINIT = 2314, 

// DCB=(RECFM=FB,LRECL=80,BLKSI ZE=-800),DI SP=NEW,KEEP), 
// SPACE=(800, (3, 1),RLSE) 

// 



Once this program has been correctly read in to the 
user's disk, it must be submitted to the operating system 
using the standard SUBMIT command: 

SUBMIT ASDF JCL PI XXXX 

XXXX is the user's computer center number. If the SUBMIT 
function Is working the following messages will be returned: 

** CARDS XFERED TO 0098P ** 

FNAME FTYPE PI SENT FOR OS BATCH JOB 
After a short delay the system will send a second message: 
FROM 0098P33 : JOB SENT TO OS WITH 21 CARDS VJRITTEN (Time) 
The user should then revert back to section B above which 



^A deck containing this program Is available from 
Professor Kirk, Sp 302. 
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discusses the TOCP command. After giving the system a period 
of time to execute the JCL job which ivas just submitted, 
execution of the TOCP command will indicate that the ASDF 
files are present. If the operating system has not completed 
executing the transfer, the TOCP request v/ill respond as in 
paragraph two above. 

B. LOADING ASDF 

Once ASDF resides on disk OSXOOl, the remaining required 
steps are simple. The user issues the command: 

EXEC ASDF191 

The system will type a series of lines onto the terminal. 
At the conclusion of these typewritten lines the user is 
prompted to issue the next two commands: 

LOGIN 192 P 

The system responds: 

192 REPLACES P (191) 

Tlie next required entry is: 

EXEC ASDF192 

These two commands v/ill cause a very long list of CP/CflS 
commands to be typed on the screen. Each of these CP/CMS 
commands is automatically executed by the operating system 
in the process of loading the ASDF programs. There are no 
furthrer commands required from the user to' complete the 
loading process. 
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'When the entire ASDF program package Is loaded, the user 
v/lll be presented with a title page and the instruction: 
"HIT SPACE AND RETURN TO CONTINUE." Once the user has typed 



"space" and "return" the display will write out 
the commands at his disposal. These commands 
been discussed In chapter four. 



a 1 I s 1 1 ng of 
have a 1 ready 
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VI . USI MG ASDF 



Each of the following sections constitutes an example 
for one of the three input modes. Each example uses only 
one input mode. In general, these modes may be combined 
during the generation of any particular filter. The mode 
selected should be the one most convenient to the user. 

A. EXAMPLE ONE - INTERACTIVE GRAPHICS 

This first example uses the ASDF program in the 

interactive graphics mode. The objective is to observe the 
responses of a digital filter which has two real zeros, one 
at the origin and another at +0.8, two complex poles located 
at (0.8 +/- jO.2), and the filter gain is 1.000. The 

interactive graphics mode is used since the exact locations 
of the poles and zeros are not critical. For the purposes of 
this example, it is assumed that the user has already 
followed the login procedure described in the above 
paragraphs, and has loaded the ASDF program. The user has 
already seen the title page, and the two pages of command 
descriptions. The screen is clear except for the top left 
corner which shows: 

♦ASDF-RE 

This group of letters is the signature for the ASDF program 
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and Indicates that the program Is ready to execute the 
user's next command. 

Since we must first enter the poles and zeros Into the 
z-plane the appropriate command to type Is: 

POLZRO 

The ASDF program will then draw the screen as shown in 
figure 6-1, and generate an audible tone. Initially the 
cursor will be located at some arbitrary position or perhaps 
may not appear on the screen at all. Movement of the thumb 
wheels moves the cursor position. Since we v/lsh to execute 
in the Interactive graphics mode/ to add a real zero, vie 
move the cursor (Intersection of the two lines tracing 
across the screen) Into the box In the left column labeled 
ARZ/ (Add a Real Zero) figure 6-2. Vie then type any 

keyboard character except "return".^ The trace cursor will 
disappear/ and the flashing alphanumeric cursor will appear. 
The program will print a carat pointing to the box selected 
by the user. The flashing cursor will move Into the box 
selected . 

The user then types the command "space"/ "e"/ "return". 
The flashing cursor disappears/ and the trace cursor 
reappears. The absence of an audible tone at this time 
Indicates that the trace cursor Is to be positioned within 
the uni t circle. 

^The 4012 terminal in the computer center requires a 
"control-s" after every response. 
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UNIT CIRCLE COMMANDS 
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Figure 6-1 Initial Screen Representation After Issuing The Command POLZRO 
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Figure 6-2 Trace Cursor Position For Executing 

The Command "Add A Real Zero"-Interactive Graphics Mode 



The trace cursor is adjusted to the oosition of the desired 
real zero, (the origin in this case), as shown in figure 
6-3. Any keyboard character except "return” is then typed. 
The ASDF program plots the zero at the origin, and returns 
the trace cursor with the audible tone. The audible tone 
reminds the user that the trace cursor is to be positioned 
in the command array. Since another real zero i s to be 
entered, the cursor is positioned as in figure 6-lb again, 
and the same steps are followed. When the trace cursor 
reappears without the audible tone, position it within the 
unit circle as near as possible to the point (0.8 +/- jO.O) 
as shown in figure 6-4. (To enter the zero precisely at 0.8 
the rectangular or polar mode would be used.) The steps for 
entering the complex pole pair exactly parallel those for 
the previous entries, and are depicted in figures 6-5, and 
6 - 6 . 

Once all the poles and zeros have been entered, the 
trace cursor may be aligned with any command box, and a 
character typed (not "return"). Once the alphanumeric 
cursor enters that box, the user types "space", "x", 
"return". The screen is cleared and the user is asked if he 
wishes to change the default filter gain. Since the default 
value is 1.0, the user responds "n", the ASDF program clears 
and responds with: *ASDF-RE 

The next step in the normal sequence is for the user to 
look at the filter characteristics. To display these 
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UNIT CIRCLE COMMANDS 
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Figure 6-3 Trace Cursor Position Required To Place A Real Zero At The Origin 



UNIT CIRCLE COMMfiNDS 




QC 

Qd 

QC_1 

CEO 

OCl. 

CD 



□□□□□□□□ 



^UJ □□□□□□□□ 

QZ 



1 LU 
QC> 
LU— . 
H-H- 
dCJ 
*— CE 



□□□□□□□□□□ 



tn rvj Q_ rvj Q_ rvi cu rvi cu cj I— 
Q CJ CJ CC QE CJ CJ QE QE ID cn 
2 : CC(XCECECDaCDQd_l 

o 



56 



Figure 6-4 Trace Cursor Position Required To Put A Real Zero At (+0.8) 
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Required To Select The Command 
"Add A Complex Pole"-Interactive 
Graphics Mode 



UNIT CIRCLE COMMANDS 
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Figure 6-6 Trace Cursor Position Required To Enter A Complex Pole At (0.8+/-J0.2) 



character i s t i cs the command typed Is: 

RESPONSE 

To ",et a hard copy of the filter responses the command; 

EXEC HRDSCPY 

is typed. Figures 6-7a, c, d, e, f, and g show the 

pertinent portions of the output. 

3. EXAMPLE TWO - RECTANGULAR MODE 

This second filter example is taken from [Ref. 5], page 
222. The filter is a fourth-order low oass Chebyshev. The 
derivation of the transfer function is described in the 
text. The given transfer function is factored to give the 

desired pole/zero locations. The filter has two complex 

zeros at (-1.0 +/- jO.O) and two complex pole pairs, one 

pair at (.7498 +/- j.5348), and one pair at (.7774 +/- 

j.2120). The filter gain Is given as 0.001836. 

In a manner exactly paralleling the previous example, 
the user selects the command ACZRCT (which stands for Add a 
Complex Zero - Rectangular Coordinates), by positioning the 
trace cursor in the appropriate box, as in figure 6-8a. 
Once a keyboard character Is typed, and the alphanumeric 
cursor moves Into the box, the "space", "e", and "return" 
keys are typed. The ASDF program will erase the screen and 
drav/ the display shown in figure 5-8b. The user then tyoes 
in the rectangular coordinates of either complex zero in 
the first complex zero pair, figure 6-8c. The system will 
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Figure 6- 7a Pole Zero Locations For Example Two 
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FILTER POLE ZERO LOCRTIONS 
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Figure 6- 7c Unit Sample Response For Example Two 
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Figure 6-7d Step Response For Example Two 
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Figure 6-7e Phase Response For Example 
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Figure 6-7f Magnitude Of H(z) For Example 
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Figune 6-7g Magnitude Of The Transfer Function In Decibels For Exampl 
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Figure 6-8b Screen Display For The Command ACZRCT 
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Figure 6-8c Screen Display After Entering The First Complex Zero 



ask whether or not the root location has been correctly 
entered. The user responds *y' for yes, *n' for no. ‘ If the 
root has been incorrectly entered, an *n’ response will 
rewrite the screen for another attempt. Upon receiving a 
'y* response ASDF rewrites the unit circle and command 
array, and plots the new root. The second complex zero is 
entered exactly as the first. 

The complex poles are entered as were the zeros. The 
trace cursor is first placed in the box ACPRCT. (Add a 
Complex Pole - Rectangular Coordinates), as In figure 6-9. 
Figure 6-10 shows the unit circle after all poles and zeros 
have been entered. The POLZRO command is terminated by 
aligning the trace cursor with any box, typing a character, 
waiting for the flashing alphanumeric cursor to enter the 
box, and typing "space”, "x". The system responds by asking 
about the filter gain. In this example the filter gain is to 
be changed to .001836, so the correct response is 'y*. The 
new filter gain, 0.001836, is typed to terminate the POLZRO 
command . 

To view the filter responses, the user enters the 
command : 

RESPONSE 

when the *ASDF-RE signature appears. The hard copy 
responses of the filter are shown in figures 6-lla through 
6-llg. 
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Figure 6-9 Trace Cursor Position To Execute Command ACPRCT 
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Figure 6-10 Unit Circle After All Poles and Zeros Have Been Returned-Example Two 
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Figure 6- 11a Pole Zero Locations For Example Two 
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FILTER POLE ZERO LOCATIONS 

Figure 6- 11b Pole Zero Locations For Example Two 
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Figure 6-llc Unit Sample Response For Example Two 
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Figure 6-lld Unit Step Response For Example 
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Figure 6-lle Phase Response For Example Two 
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Figure 6-llf Magnitude of HCz) For Example 
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Figure 6-llg Magnitude Of H(z) In Decibels For Example 
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C. EXAMPLE THREE - POLAR MODE 

This example uses the polar mode to ;^enerate an all pass 
filter. The filter complex pole pairs are to be located at 
(0.2 at +/- 45 degrees)/ and (0.5 at +/- 30 degrees). The 
filter zeros are to be at (5.0 at +/- 45 degrees)/ and (2.0 
at +/- 50 degrees). The filter gain is to be 0.01. The 
first step is to change the specifications in degrees to 
radians. Thirty degrees equals .5235 radians. Forty five 
degrees is .7854 radians. The entry of ooles and zeros 
exactly parallels the above examples/ however/ the desired 
comands are ACPPLR/ (Add a Complex Pole - Polar 
Coordinates)/ and ACZPLR/ (Add a Complex Zero - Polar 
Coordinates). Figure 6-12a shows the trace cursor position 
for ACPPLR. Figure 6-12b shows the trace cursor position 
for ACZPLR. The angular part of any pole or zero entered in 
the polar mode must be in radians. Figure 6-12c shows the 
entry of the first complex zero pair. Figure 6-12d shows 
the entry of the first complex pole pair. Figure 6-12e 
shows the change in filter gain from 1.00 to 0.01. Figures 
6-13a through 6-13g show the hard copy responses for this 
filter. It should be noted that zeros which plot outside of 
the unit circle are indicated by stars. All such zeros are 
plotted at a radius of 1.1 along the same angular component 
as the actual zero location. Zeros are plotted this way for 
two reasons: first to make the output display consistent for 
all relevant filters/ and second to emphasize zeros v/hich 
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UNIT CIRCLE COMMANDS 
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Figure 6-12a Trace Cursor Position For The ACPPLR Command 
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Figure 6-12b Trace Cursor Position For The ACZPLR Command 



o 

o 

o 

- 3 * 

m 

00 

o 










o 


CO 






• 








• 






UJ 


X 


CO 


% 


X 


• 






O- 






c 


X 


X 




CJ 


CO 






O X 




c 






< 


< 




UJ 






o- 


UJ 




> 


c 




X 


CO 


X 






o 


X 


< 


c 




Ci 


— 




c 






o H* 


h- 




a: 


a: 


< 


CJ 


CO 


CG 






o c 






a. 






UJ 


3 




UJ 




O X 


O 


< 




UL 




c: 


X 


UJ 


_l 




• 


H 






o 


< 






X 


CJ 




in UL 




c: 






H- 


< 


X 


H- 


X 


in 


i-H — 


O 


< 


X 


UJ 


UJ 








< 


CO 




cr 


d: 




o 


X 


CO 


< 


C 




r-* 


UJ % 


UJ 




— 


X 


h- 


UJ 




H- 




• 


cc- 


M 


UJ 




H- 




cc 


UJ 








L >- 




X 








— 


h- 


O 




r 


X 


h- 


CO 


X 


X 


X 


< 


UJ 




1 


_i 


UJ 




H- 


CJ 


UJ 




— 


X 






<1^ UJ 


_1 


cr: 


X 


< 


1- 


UJ 


cc 








X CL 


o. 


UJ 


UJ 






cr 


CL 


LL 


UJ 




>- 


X 


h- 


X 




X 




c 


X 


Cl 




CO K- 


c 


X 


c 


UJ 


< 


q: 


cs: 


o 






O 


o 


UJ 


C- 


X 


X 


UJ 


a. 


CJ 


H- 


o 


CC h- 








H- 


H* 


CG 


a. 




— 


• 


UJ CJ 


< 


V 


c 








< 


UJ 


X 


in 


M UJ 




>- 


CJ 




CO 


X 




cc. 


c 


rH 


CC 


o 


< 




• 


CO 


X 


CO 


< 


< 




X C' 


c 




< 


CO 


UJ 




< 




X 




UJ o 


< 


a. 


H 


UJ 




T* 




X 




“ 


_l CJ 




</> 


UJ 


X 




CJ 


c 


H* 






o. 


c 




X 


o 


UJ 


< 


X 


c 




/N 


X CL 


h- 


a. 


H- 


cc 


CO 


UJ 


< 


CG 






c — 



CJ 



3 * 

UJ 



83 



Figure 6-12c Screen Display After Entering The First Complex 
Zero Pair In Polar Coordinates 



TO ADD A COMPLEX POLE WITHIN THE 













X 












(V 










o 


r— 










o 


o. 










o 


E 










-a* 


c 










LTk 


o 










OC 














4-> 










• 












o 


L- 










¥ 


LL 






















+ 


c 


tu c o 


• CO • 






• 


X 


X 2 tu 


CO X sX • 






o - 


1- 


1— < c 


0^ < -J CD CO 






O X 




— 


UJ — <: — UJ 






o- 


b£ 


z -J > 


-J C X X 






o 


C 


— < c 


< — c 






o 1- 




X — c:: 


ct: a: o (/:• cc 


UJ 




o c 




1- c a. 


O UJ 




j- 


O X 




— < 


X o X UJ 


CD 


lA 


CVJ 


4-» 


X 


UJ — — X 


X 


oo 


• u. 


c 




X < X »- 


< 




o — 


UJ 


tu UJ X 


O C2 




• 






-1 X »- 


UJ CO < c 






UJ X 


u 


O 1 


UJ q: UJ 1— 






CD - 




CL S 


CO UJ ci UJ 






L >- 




a: 


1 1-0 




^ -J- 


4- 


X UJ CO 


1- X X < UJ 




1 


-J 


< 


UJ K 1- 


00 UJ c' — X 






UJ 




-J X X 


X ujcr: — 




+ 


X a. 


> 


Q. UJ UJ 


X UJ a: a. LL 


UJ 




>• 


03 


X 


CQ O ^ 


o 




00 »— 




o UJ c 


O cr cr: o 


X 




UJ 


C. 


o -J a. 


X K UJ CL o 


h- 




— J h- 


{/) 


O X 


cr CO CD a. 


— 


CNJ 


c o 


•— 


< a: o 


X X < UJ 


X 


• 


a. UJ 


o 


— o 


XX a: 


a 


c 


a: 




o o 


• X CO < 


< 




X o 


c 


c: < 


00 < < 


X 




UJ o 


a; 


< 1- K 


UJ 1- X X 






-J o 


(U 


— UJ 


XUJOO»— 






a. 


u 


O X X 


C X < X o 






X u- 


u 


I- X H 


CD 1— UJ < CD 






o — 
o 


00 












X 










UJ 












X 


r~* 

1 

CC 










UJ 










X 












K 


o 












u 












3 












bC 












•— 












LL 



84 



Pole Pair In Polar Coordinates 



THE FILTER GAIN IS CURRENTLY 1.00000000 



c 







00 




















LU 










CQ 










X 






0- 




X 










X 




o 


— 








o 


< 




cr: 




o 


c: 




LU 




o 






3* 




o 


c: 




LU 




o 


LU 




u. H 




rH 










O 






q: — 






— 




o o 




o 


ll 




O- 










CN 






tu 




l-H -J 




• • 


X 


X 


< 




i/) 


H 




X 2 








a: 


•• — 






LU 


c 


< o 






CJ 




O LU 




- 


X 


>- 


X 




< 


< 




cr 




O 


X 


LU 


LU < 


rH 




o 


Q. 




O 


X 




>- 


_J LU 


• 


LU 


c 




— Q 


O 


K 


K 




U- X 




-J 












X 




X o 




LL 






LU X 










X — 














LU 






LU 




X 


X 




X 






o 




K 






>- 




X 






c 















LU 










a. 










>- 

















85 



Figure 6-12e Screen Representation After Changing Filter Gal 
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UNIT CIRCLE 
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FILTER POLE ZERO LOCATIONS 
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Figure 6-13c Unit Sample Response For Example Three 



09 •] 



& 

0 

0 

0 

0 

0 

O- 



o 

o 



•o 

•-0 

o 

o 

oO 

« 

u=r !l 



ex 

0(J3 

O 

.ever 

— UJ 



OLl. 

o 

« 

Uo 



o- 

0 




0 



• 0 



o 

o 

CD 



0 



© 



o 

o 

'n* 



© 



0 



o 

o 

'csj 



00*1 



-0 



08*'0 0 fi *0 

3SNDdS3y d3iS 



— 0- 
00*0 



OTi*0 



o 

o 

• 

o 



89 



/ 



Figure 6-13d Unit Step Response For Example Three 
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Figure 6-13e Phase Response For Example Three 
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Figure 6-13f Magnitude Of H(z) For Example Three 
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Figure 6-13g Magnitude Of H(z) In Decibels For Example Three 



are close to, but outside, the unit circle. The importance 
of distinguishing between zeros very near the unit circle. 



either interior or exterior, rests on the 
sensitivity of the phase of the frequency response as 
moves from the interior, onto, and to the exterior 
unit circle. 



extreme 
a zero 
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VII. ADDITIONAL EXAMPLES 



The purpose of this chapter is threefold. The primary 
objective is to provide the user with a series of plots 
demons t rat i the general form of those types of filters not 
considered in chapter six. Additionally, this chapter 
demonstrates the responses /generated by poles located at a 
series of positions within the unit circle. It further 
provides examples for user comparison when learning to 
interact with the z-plane. 

The filters considered in the examples of chapter six 
include a fourth-order low-pass filter and a second-order 
all-pass filter. The sequence of figures 7-la through 7-lg 
shows the general pole zero placements for a high-pass 
filter, and associated responses. Similarly, fitrure 7-2a 
through 7-2g show the responses for a band-pass filter. 
Figures 7-3a through 7-3g show a band-reject or notch 
f i 1 ter . 

The secondary purpose of the chapter is to show the user 
the general forms of time and frequency responses as the 
filter poles move in the z-plane. The sequence of figures 
7-4 through 7-11 shows the responses for a real pole at 
various locations on the real axis. Figures 7-12 through 
7-19 similarly show the responses of a complex oole pair 
moving across the z-plane. Frequency responses are not 
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shown for poles on, or otitside of, the unit circle since the 
frequency response is not defined for these cases. Figures 
7-17, 7-20, and 7-21 shov/ the change in responses for a pair 
of complex poles which move away from the real axis. 

The final intent of the chapter is to provide the user 
with additional figures for comparison should he need 
further practlve in manipulating poles and zeros in the 
z-plane. It is important to note that the ASDF program will 
not allow the user to position poles on or outside the unit 
circle. This limitation results in several plots which can 
not be reproduced by the user. 
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Figure 7-la Pole Zero Locations For A Second-Order High-Pass Filter 
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FILTER POLE ZERO LOCATIONS 
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Figure 7-lc Unit Sample Response For A Second-Order High-Pass Filter 
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Figure 7-ld Unit Step Response For A Second-Order High-Pass Filter 
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Figure 7-le Phase Response For A Second-Order High-Pass Filter 
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Figure 7-lf Magnitude Of H(z) For A Second-Order High-Pass Filter 
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Figure 7-lg Magnitude Of H(z) In Decibels For A Second-Order High-Pass Filter 
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Figure 7-2a Pole Zero Locations For A Fourth-Order Band-Pass Filter 
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Figure 7-2c Unit Sample Response For A Fourth-Order Band-Pass Filter 
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Figure 7-2d Unit Step Response For A Fourth-Order Band-Pass Filter 
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Figure 7-2e Phase Response For A Fourth-Order Band-Pass Filter 
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Figure 7-2f Magnitude Of H(z) For a Fourth-Order Band-Pass Filter 
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Figure 7-2g Magnitude Of H(z) For A Fourth-Order Band-Pass Filt 
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Figure 7-3a Pole Zero ' Locations For A Fourth-Order Notch Filter 
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FILTER POLE ZERO LOCATIONS 
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Figure 7- 3c Unit Sample Response For A Fourth-Order Notch Filter 
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Figure 7- 3d Unit Step Response For A Fourth-Order Notch Filter 
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Figure 7-3e Phase Response For a Fourth-Order Notch Filter 
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Figure 7-3f Magnitude Of H(z) For A Fourth-Order Notch Filter 
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Figure 7-3g Magnitude In Decibels For A Fourth-Order, Notch Filter 
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FILTER POLE ZERO LOCATIONS 
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Figure 7-4b Unit Sample Response For A Single Real Pole At (-1.05) 
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Figure 7-4c Unit Step Response For A Single Real Pole At (-1.05) 
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Figure 7-5b Unit Sample Response For A Single Real Pole At 
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Figure 7-5c Unit Step Response For A Single Real Pole At 
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Figure 7-6b Unit Sample Response For A Single Real Pole at 
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Figure 7-6c Unit Step Response For A Single Real Pole At 
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Figure 7-7d Phase Response For A Single Real Pole At 
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Figure 7-6e Magnitude Of H(z) For A Single Real Pole At 
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Figure 7-6f Magnitude Of H(z) In Decibels For a Single Real Pole at 
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Figure 7-7b Unit Sample Response For A Single Real Pole At 
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Figure 7-7c Unit Step Response For A Single Real Pole At 



o 

ojo 




132 



Figure 7-7d Phase Response For A Single Real Pole At 
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Figure 7-7e Magnitude Of H(z) For A Single Real Pole At 
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Figure 7-7f Magnitude Of H(z) For A Single Real Pole At 
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FILTER POLE ZERO LOCfiTIONS 
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Figure 7-8b Unit Sample Response For A Single Real Pole At (0.3) 
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Figure 7-8c Unit Step Response For A Single Real Pole At (0.3) 
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Figure 7-8d Phase Response For A Single Real Pole At (0.3) 
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Figure 7-8e Magnitude of H(z) For A Single Real Pole- At (0.3) 
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Figure 7-8f Magnitude Of H(z) In Decibels For A Single. Real Pole At (0.3) 
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FILTER POLE ZERO LOCATIONS 



•I 

o 

a 

It 



o CL 

Gio^ 

(L^cn 

^CMUJ 



Q ^ 
ou. 
o 

Gh*^ 

^CM 



Q 

G- 

G- 

G 



G— 



G- 

G— 



o 

o 



.CD . 



o 

o 

• 

,c\J 



O 

o 

'oo 



o 

o 



o 

o 



00‘I 08*0 09*0 OTi'O 02*0 OO^cP 

3SN0dS3y 3Sindwr 



14 2 



Figure 7-9b Unit Sample Response For A Single Real Pole At (0.8) 
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Figure 7-9c Unit Step Response For A Single Real Pole At (0.8) 
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Figure 7-9d Phase Response For A Single Real Pole At (0.8) 
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Figure 7-9e Magnitude Of H(z) For A Single Real Pole At (0.8) 
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Figure 7-9f Magnitude Of H(z) For A Single Real Pole At (0.8) 
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FILTER POLE ZERO LOCRTIONS 

Figure 7-lOa Example Of A Single Real Pole On the Unit Circle at (1.0) 
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Figure 7-lOb Unit Sample Response For A Single Real Pole At (1.0) 
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Figure 7-lOc Unit Step Response For A Single Real Pole At ('l.O) 
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FILTER POLE ZERO LOCATIONS 
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Figure 7-llb Unit Sample Response For A Single Real Pole At (1.05) 
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Figure 7-llc Unit Step Response For A Single Real Pole At (1.05) 
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FILTER POLE ZERO LOCATIONS 
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Figure 7-12b Unit Sample Response For A Complex Pole Pair At (-0.8 j;^j0.5) 
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Figure 7-12c Unit Step Response For A Complex Pole Pair At (0. 
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FILTER POLE ZERO LOCATIONS 
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Figure 7-13b Unit Sample Response Of A Complex Pole Pair At (-.866025 +j0.5) 
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Figure 7-13c Unit Step Response Of A Complex Pole Pair At .(-0.866025 j^jO.5) 
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FILTER POLE ZERO LOCRTIONS 



a 

o 



‘ OtT0 



09’1 08‘0 

3SN0dS3y 



OQ 

a 

o 

o 

o . 



-3* tl 



aCC 

0(_3 

o’ 

.ojQC 

-•UJ 



OLl. 

o 

.a 



00-0 08 - 0 - 

3sindwi 



o 

o 



.a . 

QO - 



O 

O 

• 

.a 

CD 



O 

O 

.a 

=f* 



o 

o 

• 

.a 

OJ 



o 

o 



lO‘ 

m 

o 

•m 

+ l 

00 

o 

\ 



09-1- 






160 



igure 7-14b Unit Sample Response For A Complex Pole Pair At 
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Figure 7.14c Unit Step Response For a Complex Pole Pair At 



o 



OJQ 




16 2 



Figure 7-14d Phase Response For A Complex Pole Pair At (-0.8 +j0. 
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Figure 7-14e Magnitude Of H(z) For A Complex Pole Pair At 
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Figure 7-14f Magnitude Of H(z) In Decibels For A Complex Pole Pair At (-0.8 +j0.5) 
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FILTER POLE 'ZERO LOCFUTONS 

Figure 7-15a Example Of A Complex Pole Pair At (-0.3 +j0.5 
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Figure 7-15b Unit Sample Response For A Complex Pole Pair At 



OOM 08*0 09-0 0fi‘0 

3SN0dS3y d3iS 



02-0 



o 

o 

•o 

'"o 

o 

o 

• 



CE 

OCJ3 

.cuce: 
— LlJ 



OLl. 

o 



o 

o 

'oo*: 



o 

o 

• 

CD 



O 

O 



O 

O 



O 

O 



00 ‘0 



LO 

o 

•m 

+ 1 

CO 

o 

I 



167 



Figure 7-15c Unit Step Response For A Complex Pole Pair At 
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Figure 7-15d Phase Response For A Complex Pole Pair At 
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Figure 7-15e Magnitude Of H(z) For A Complex Pole Pair At 
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Figure 7-15f Magnitude Of (H:^ In Decibels For A Complex Pole Pair At (-0.3 +j0.5) 
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FILTER POLE ZERO LOCRTIONS 
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Figure 7-16b Unit Sample Response For A Complex Pole Pair At (0. 
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Figure 7-16c Unit Step Response For A Complex Pole Pair At 
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Figure 7-16d Phase Response For A Complex Pole Pair At (0.3 ;tjo.5) 
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Figure 7-16e Magnitude Of H(z) For A Complex Pole Pair At (0.3 +j0.5) 
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Figure 7-16f Magnitude Of H(z) In Decibels For A Complex Pole 
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FILTER POLE ZERO LOCATIONS 
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Figure 7-17b Unit Sample Response For A Complex Pole 
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Figure 7-17c Unit Step Response For A Complex Pole Pair At (0. 
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Figure 7-17d Phase Response For A Complex Pole Pair At (0. 
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Figure 7-17e Magnitude Of H(z) For A Complex Pole Pair At (0. 
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Figure 7-17f Magnitude Of H(z) In Decibels For A Complex Pole 
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FILTER POLE ZERO LOCATIONS 

Figure 7.18a Example Of A Complex Pole Pair On The Unit Circle At (0 . 8S6025*i0 . 5) 
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Figure 7-18b Unit Sample Response For A Complex Pole Pair At (0.866025 +j0.5) 
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Figure 7-18c Unit Step Response For A Complex Pole Pair At (0.866025 +j0.5) 
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FILTER POLE ZERO LOCUTIONS 
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Figure 7-19b Unit Sample Response For A Complex Pole Pair At (0.9 +j0.5) 



o 

o 




188 



Figure 7-19c Unit Step Response For A Complex Pole Pair At (0.9 +j0.5) 
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FILTER POLE ZERO LOCATIONS 
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Figure 7- 20b Unit Sample Response Of A Complex Pole Pair At (0.8 +_j0.1) 
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Figure 7-20c Unit Step Response For A Complex Pole Pair At (0.8 j^jO.l) 
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Figure 7-20d Phase Response For A Complex Pole Pair At (0.8 +j0.1) 
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Figure 7-20e Magnitude Of H(z) For A Complex Pole Pair At (0. 



o 



OJQ 




194 



Figure 7-20f Magnitude Of H(z) In Decibels For A Complex Pole Pair At (0.8 +^j0.1) 
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FILTER POLE ZERO LOCATIONS 
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Figure 7-21b Unit Sample Response For A Complex Pole Pair. At (0. 
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Figure 7-21c Unit Step Response For A Complex Pole 
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Figure 7-21d Phase Response For A Complex Pole Pair At (0. 
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Figure 7-21e Magnitude Of H(z) For A Complex Pole Pair At(0. 
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Figure 7-21f Magnitude Of H(z) In Decibels For a Complex Pole 



VIII. INTERPRETATION OF ASDF OUTPUT 



When 


the user 


initially confronts 


the ASDF 


program i t 


is expected that 


he 


will attempt to duplicate 


some 


of the 


exampl es 


presented 


earlier in this 


thes i s 


or 


s i mi 1 a r 


exampl es 


found 


t n 


the 1 i terature . 


These 


types of 


appl i cat ions of 


ASDF are straightforward and 


the 


resul ts 


generated 


can 


be 


compared wi th 


pub 1 i shed 


resul ts . 



Interpretation of the output plots is seldom confusing. 

As the user becomes more familiar with the use of the 
ASDF program and the positioning of the poles and zeros in 
the z-plane he will undoubtedly attempt to generate filters 
with sharper transition regions, and greater attenuation in 
the stop band(s). In general these filter improvements are 
accomplished by increasing the order of the system and 
positioning the poles closer to the unit circle. As the 
order of the system increases, and as the poles of the 
filter become closer to the unit circle, the filter becomes 
increasingly subject to the difficulties associated with 
the finite precision arithmetic used in the filter 
impl ementat ion . 

A tacit assumption made very early in most digital 
signal processing courses (and often quickly forgotten) is 
that infinite precision arithmetic is available in 
performing the calculations. This assumption breaks down 
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when digital filters are realized in hardware or simulated! 
In general there are three sources of error which may affect 
the performance of digital filters. Due to the finite 

precision representat ion of the coefficients, the values of 
the implemented filter coefficients are only approximately 
equal to the desired values. Second, errors may be generated 
in the arithmetic operations. In general, the product of 
two n-bit numbers results in a 2n-bit result. An error 
occurs when this 2n-bit number is converted to n-bits for 
succeeding computations. The third type of error 

encountered is the analog to digital conversion required for 
processing analog signals with digital filters. Only the 
first two of these errors is of concern in the ASDF program. 

Vlhile an in depth analysis of finite precision effects 
is beyond the scope of this work, let us consider the filter 
v/hich has eight zeros located at the origin, and eight poles 
located at -0.9. Entering these locations into the z-plane 
using the ASDF program generates the filter responses in 
figures 6-la, b, c. A superficial glance at these responses 
Indicates that this filter is unstable. The filter must be 
stable, however, since all of the poles lie within the unit 
circle. 

There are several phenomena involved in this example. 
The first is the dynamic range limitation of the system on 

which the filter is implemented. The IBM-360 is limited to 

6 U’ 6 3 

computing with numbers between 15 and -16 . Once the 
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FILTER POLE ZERO LOCATIONS 
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Dynamic Range Limitations 
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Figure 8-lc Step Response For A Filter Exceeding System Dynamic Range Limitations 



larger of these limits is exceeded, the computations are 
terminated and it is not clear whether the designed filter 
is actually unstable or has merely exceeded the dynamic 
range limitations of the machine. The dedicated user may 
decide to modify the programs presented in the appendicies 
to investigate the time responses of such systems. Clearly 
the form of the time responses will not be significantly 
altered if a modified unit sample response 

1 X 0, 0, 0, . . . 

and a modified step response 

. in-10 1 irt-10 1 in-lO 

1 X 10 , 1 X 10 , 1 X 10 , . . . 

is used. 

If these modifications are performed, two other 

phenomena become apparent, each relating to the finite 
precision implementation of the digital filter. The first is 
that, in general, digital filters have deadbands, the second 
that the filter output may degenerate to a limit cycle. The 
original works by Blackman [Ref. 1] . demonstrate that 
digital filters exhibit both of these effects. Extensive 
covgerage of these phenomena is not appropriate here, 
however, a dead-band in a digital filter can be 

s impl i s t i ca 1 1 y described as a span of filter input values 
for which the filter output does not change. The limit 
cycle effect is evident when a digital filter responds to a 
steady state input. For constant inputs, the filter output 
may run through a cycle of values. The propensity for 
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exper I enc i nj; limit cycles is obvious in the calculation of 
the step response. A close look at the unit sample response 
shows that after a sufficiently lar<te number of input 
values, this unit sample sequence becomes in reachine: 
steady state, a constant. 

When observing the output response of the example 
filter, there is still another possibility for explaining 
this output. This explanation revolves around the fact that 
the poles in the implemented filter are not at the locations 
expected. In fact, the poles may actually move outside the 
unit circle even though they were entered within. In 
general, there are only a finite number of pole locations 
possible with finite word length coefficient 
representat i ons . The effect of this limitation is best 
explained using another example. Consider a filter with 
three zeros located at 0.5, and 0.2 +/- JO. 5, and four poles 
located at -0.999, -0.6, and -0.69 +/- JO. 69. The infinite 
precision transfer function for this filter is expressed in 
equation 8.1. 

- .9 + .49 Z - .145 (8.1) 

Z*^ + 2.979 Z^ ♦ 3.75822 t } + 2.3497398 Z + .57074868 

Now consider the finite precision case where each of the 
coefficients must be represented by a six bit number. The 
pole at -0.6 is expressed in binary (exclusive of sign) as 



207 



.10011001100. . . and Is truncated to six bits as .100110. 
i'ihen all of the coefficients are so represented, the 
transfer function is expressed by equation 8.2. Clearly the 
poles and zeros of this second transfer function are not the 

- .890625 +.474275 Z - .140625 (8.2) 

:(z) = ^ ^ 

Z + 2.9375 Z + 3.75 Z ♦ 2.3125 Z + .5625 

sane in equation 8.1. In j^eneral the differences will change 
all of the filter responses. In fact for this particular 
case the nev/ poles of the system are at -0.70522 +/- 
JO. 16677, and -0.76353 +/- JO. 698616; this filter is 

uns tab 1 e ! 

The reason for the shift in pole locations is clearly 
portrayed by the filter whose character 1 s t ic equation is 
equation 8.3. 

(Z + .995) ^ = 0 (8.3) 

This filter has eight poles at -0.995. Note what happens to 
the filter when an error is Introduced by truncating the 
infinite precision coefficient values to thirty-two bits. 

For simplicity restrict the discussion to the ramifications 

0 

of an error generated in the value of the constant (.995) 
This error can be as large as 2 x 10 With this error 

explicitly included, the new characteristic equation will 
be : 

(Z + .995)® + 2 X 10”^°= 0 (8.4) 
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The roots of equation (3.4) are complex, and 1 1 e on a circle 
of radius 0.06 around the point -0.995 +/- J 0.0. Figure 
8-2 shows that three of these roots lie outside the unit 
circle. Thus with a finite precision implementation of the 
filter (even v/i th 32 bits of precision), it is possible for 
the pole locations to move so far that the filter actually 
becomes unstable. 

The results of this section as well as four types of 
ASDF outputs which have been observed, are summarized in 
figure 8-3. Trace A shows the unit sampie response of a 
stable filter. Trace B shows the response of a filter which 
exceeds the dynamic range limitations of the implementation 
system. Trace B1 is for an unstable filter and B2 is for a 
stable filter caught in a limit cycle. These two traces can 
be observed or explored only by modifying the unit sample 
sequence. Trace C shows a stable filter which falls into a 
limit cycle below the dynamic range limit of the machine 
upon which the implementation has been performed. Note that 
the performance of a filter above the dynamic range 
limitation of the IBM 369 (or any other machine) cannot be 
observed, unless the ASDF program is modified. 
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Unit Circle 
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Figure 8-2 Three Of These Complex Poles Are Outside The Unit Circle 
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Figure 8-3 Sununary Of Possible Observed Unit Sample Responses 
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IX. SUMMARY AND CONCLUSIONS 



The Advanced Simulation of Digital Filters has been 
implemented on the IBM 360, and is intended for the use by 
persons beginning their study of digital signal processing. 
The programs included within the package have been designed 
to make the user's interaction with the z-plane as simple as 
possible. Vihile the programs themselves are available for 
study, all of the computations are invisible to the user, in 
an attempt to develop the user's intuition of the discrete 
domain without overwhelming him with arithmetic 
computat i ons . 

The ASDF package provides the user with an interactive 
method by which the filter pole and zero locations can be 
manipulated and the effects observed of such movement. ASDF 
can simulate filters of up to tenth order, and is punctuated 
with numerous comments and instructions which are to assist 
the user. The interface v/i th the z-plane is available in a 
rapid interactive graphics mode, as well as in the slower 
polar and rectangular modes, each providing greater position 
location accuracy. The filter responses generated by ASDF 
are presented as a series of displays directly on the 
Tektronix 4012 terminal. Versatec plot capability has been 
integrated into the system to orovide high quality copies of 
the system responses. 
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A. CONCLUSIONS 



There are several points which are v/orth considering for 
persons Intending to pursue this or other areas of 
interactive systems development with the graphics terminals. 
While comments may apoear obvious they are stated to 
minimize the difficulties v/hich can be encountered. First, a 
serious attempt should be made to define the entire data 
structure before beginning the programming stage of the 
work. All phases of the effort must be planned before 
commencing with the actual programming to insure that the 
data structure selected is sufficiently general, and 
sufficiently flexible to endure the additions and changes 
which inevitably occur when initial objectives are 
overlooked or changed. Second, it is imperative to master 
all of the software tools available and applicable to the 
project. There are features of the PLOT 10 software which 
exist and were not used in the ASDF program. This 
shortcoming is due to insufficient a priori knowledge. For 
example, the writing of text portions of the screen 
presentations should have been much simpler. Finally, it is 
important to research all of the attributes of the devices 
used for the final problem solution, and direct the solution 
to take advantage of all the system capabilities. 
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B. AREAS FOR FUTURE WORK 

There are several areas in which this orogram could be 
improved or extended. The possibility exists for 
transferring the ASDF programs to the Tektronix 4082 
computing system. This move would avoid many of the 
shor tcommi ngs associated with the current IBM 360 system. 
The user would be relatively immune to system crashes, and 
the constraints placed upon him by concerning the IBM 
system's being operational and/or available. There are 
several disadvantages associated with this move. With the 
advent of the new computer system the majority of the 
current problems concerning system reliability, and 
constraints on available access hours should be eliminated. 
Furthermore, transferring these programs to the 4082 system 
will require significant rewriting to make execution 
possible on a system which has a much smaller virtual 
machine. Furthermore to fully benefit from the 4082 system, 
the Interfacing routines should be rewritten to take 
advantage of the 4082 refresh capability. It seems that all 
of the effort involved for this change is not worth the 
expected improvements. 

The ASDF program is now useful for analyzing digital 
filters. It seems likely that this capability could be 
expanded into a method for computer-aided design of digital 
filters. There are standard transformations which allow the 
filter designer to mao analog filter requirements into the 
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appropriate digital requirements, and as well there are 
transformations which make possible the computation of 
digital filter transfer functions from known analog transfer 
functions. There are also a variety of optimization 
techniques referenced in the literature which could be 
i mplemented. 

The other area of obvious expansion for the ASDF program 
is to include the capability for oerforming spectral 
analysis. It would be worthwhile (though perhaps not 
particularly Involved) to include in the system the ability 
to calculate Discrete Fourier Transforms, and graphically 
display the results. Additionally, it would be possible to 
demonstrate the correspondence betv/een performing 
convolution in the time domain directly and performing the 
convolution in the frequency domain by multiplying the 
DFT's. The relation between the circular and linear 
convolution could also be graphically portrayed. 

As a final precaut ionary note, it is not possible to 
overestimate the length of time required to generate well 
documented, correct, useful programs on an interactive 
basis. Changes in system parameters, and the limitations in 
system softv/are can make the process of simplifying the 
user's interactive responsibilities a very difficult process 
for anyone not thoroughly acquainted with all available 
resources. The myriad of possible incorrect user responses 
which must be anticipated make the problem extensive. 
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APPENDIX A: PROBABLE SOLUTIONS FOR DIFFICULTIES 



SYMPTOM: RESPONSE command returns an error message upon 

execution. 

SOLUTION: Check the virtual machine size with the command 
"CP Q F" to insure a virtual machine size of LOOK. If the 
user has only been allocated 256K, logoff and log back on as 
directed in chapter five. Load ASDF as described in chapter 
six. 

SYMPTOM: User has inadvertently entered the control 

program, ie., CP. 

SOLUTION: Type in "b" and "return." Normal program 

execution should resume. 

SYMPTOM: The carat appears with the blinking cursor, but it 

is not located in the command array. 

SOLUTION: Move the trace cursor into the command array, and 

type a "space." Normal operation should resume. 

SYMPTOM: No terminal activity. 

SOLUTION: First, check to see whether or not the system has 

crashed. If the system has crashed, logon as in chapter 

five, and load ASDF as in chapter six. If the system has not 
crashed, move the cursor to the center of the screen and 
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type "space” and "return." If normal operation has not 
resumed hit the "break" key. When the CP prompt appears, 

type "b." Normal operation should resume. 

SYMPTOM: An error Is detected In the execution of ASDF191 

EXEC Indicating that no temporary disk space Is available. 
SOLUTION: Check with the consultant or a systems 

programmer . 

SYMPTOM: The user can not escape from a delete pole or 

delete zero command In the Interactive graphics mode (lAG). 
SOLUTION: Implicit In each request to delete from the 

z-plane Is a request, validating the change to be made. In 
the delete lAG mode the trace cursor returns after 

Illustrating the root to be deleted. Type "y" If the correct 
root was Identified. Type "n" If the wrong root was 

Identified, and attempt the procedure again. 

SYMPTOM: The user did not enter the correct amount of data 

and the program terminated on a 217 error. 

SOLUTION: Type "b" and "return." The program should 

restart with the ASDF signature. At this point reenter the 
command desired. If executing the POLZRO command It Is 
advisable to execute the list (LST) command to check on the 
current filter status. 
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APPENDIX B: IMPORTANT SUBROUTINES AND VARIABLES 



The subroutines and variables described in this appendix 
are applicable to all programs within the ASDF package. 

A. SUBROUTINE NAMES 

The subroutines within the ASDF package are divided into 
three major cataj^ories: 1) Plot-10 software, 2) Versatec 

software, and 3) routines generated explicitely for the ASDF 
prof^rams . 

1 . Plot-10 Software Subroutines 

The Plot-10 subroutines are all adaquately described 
in reference 10, and are listed below. 



ERASE 


MOVABS 


ANMODE 


RECOVR 


HOME 


LINEF 


NEHLl N 


SWI NDO 


V'/INDO 


POLTRN 


MOVEA 


DRAV/SA 


LI MTRN 


DRAWR 


DRVIREL 


MOVREL 


VCURSR 


BELL 


PO 1 NT A 


SCURSR 


INIT 


TABHOR 


TABVER 


SEELOC 


FIN 


F 1 N 1 TT 


DR AW A 


DR'/ABS 
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2 . 



Versatec Subroutines 



The Versatec subroutines are listed belmv, and are 
explained in detail in reference 7. 



PLOTS 


PLOT 


Wl NDOW 


SCALE 


AXIS 


SYMBOL 


LI NE 







3 . ASDF Subroutines 

Each of the following subroutines have been 
generated for the ASDF programs. 

CRRMMM - Add and delete poles and zeros from the z-plane. 
The change "C" can be either "A" for add or "D" for delete. 
The roots "RR" can be "RZ" for real zero, "CZ" for comolex 
zero, "RP" for real pole, or "CP" for complex pole. Mode 
"MMM" can be "I AG" for interactive graphics, "POL" for 
poles, or "RCT" for rectangular. 

FLTR2 - Prints a listing of the pole and zero locations of 
the filter being considered. 

ASUBK “ Generates the aj^ coefficients required to evaluate 
the transfer function. 

3SUBR - Generates the bp coefficients required to evaluate 
the transfer function. 

ADJUST - Insures that the correct number of zeros located at 
the origin are included in the coefficients a^^ and b^. 

RSPNSE - Computes the time responses of the filter. 

FRQNCY - Computes the phase and magnitude of the transfer 
function for the prescribed filter. 
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I. 






PRPPLT 



- Scales the ordinate data prior to plottins' the 
responses . 

PLTIMP - Plots the unit sample response. 

PLTSTP - Plots the unit step response. 

PLTTRF - Plots the phase of the transfer function as well as 
the magnitude in linear and decibel format. 

UNTCR - Plots the unit circle with the appropriate poles and 
zeros . 

START - Initializes all arrays and required variables. 

FINISH - Checks to insure that the constructed filter is 
causal/ and generates a file to save the pole and zero 
locations. 

ANGLE - Draws the angle sign for polar commands. 

ERROR - Controls orogram flow and error messages for user 
generated errors. 

LST - Prints a listing of the current pole zero location. 
PLTTBL - Plots the poles and zeros on the unit circle. 

CHCKSZ - Insures that roots are available to be deleted and 
that the maximum system order will not be exceeded. 

UPDATE _ Makes changes in the POLZRO array to corresoond to 
user changes in the z-plane. 

LOCATE - Searches the POLZRO table for the root closest to 
that which the user wishes to delete. 

STRRTS - Puts poles and zeros into the POLZRO array. 

PLOTRT - Plots a single pole or zero on the unit circle. 
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B. IMPORTANT VARIABLES 

The variables described below pertain to all of the 
subroutines and pro;trans within the ASDF packafte. 

IROOTS(I) - Keeps track of the number of poles and zeros In 
the system under consideration. IROOTS(l) Is the number of 
real zeros. IR00TS(2) Is the number of complex zero pairs. 
IR00TS(3) Is the number of real poles. IR00TS(4) Is the 
number of complex pole pairs. 

POLZRO(I^J) - The locations of the zeros are stored for 
values of I from one to ten. The pole locations are stored 
for values of I from 11 throuith 20. The real part of the 
rectan^^ular representat Ion of a root location Is stored with 
J equal to one. The Imaginary part of the rectangular 
representlon for a root location Is stored with J equal to 
two. The magnitude and angle of the root's polar 
representation are stored for J values of three and four, 
respectively. P0LZR0(I,5) holds the root type. 
lERROR - Keeps track of the generation of an error. Once 
the error has been Identified the value Is set to one. When 
the error has been explained to the user the value Is reset 
to zero. 

A( I ) and B(l) - Store the values of the a^^ and b^ 
coef f 1 c i en ts . 

AA( I ) and BB(I) - Store the double precision values of the 
aj^ and b^ coefficients. 
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YI ( I ) ■ Stores the values of the unit sample response to be 
plotted. 

YS(I) - Stores the values of the unit step response to be 
plotted . 

YP( I ) - Stores the values of the phase of the transfer 

function to be plotted. 

YM(I) - Stores the magnitudes of the transfer function to be 
plotted . 

YMOB(I) - Stores the transfer function magnitudes in 
decibels to be plotted. 

EN(I) - Stores the time values against which the dependent 
variables are plotted. 

lYIPTS - Holds the number of unit sample points to be 
plotted. 

lYSPTS - Holds the number of unit step response points to be 
plotted . 

NDEN - Is the order of the denominator plus one. 

NNUM - Is the order of the numerator plus one. 

NORDER - Is the order of the system under cons i de ra t ion . 
This can be either of the numberator or denominator. 

LAST - Holds the number of points to be plotted. 

FLTRGN - Is the value that the user assigns to the filter 
gain. 

ZEROS(I) - Holds the complex coefficients which will be made 

real to generate the b coefficients. 

r 
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POLES(I) - Holds the comolex coefficients which will be made 
real to generate the a^ coefficients. 

The remainder of the variables are either obvious, or 
used merely for control of the program flow, ie. do loop 
i ndi cies, etc. 
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APPENDIX C. STRUCTURE FOR ASDF TEXT FILE 



The file ASDF TEXT is a composite file made up of all of 
the required programs and EXEC files for running the ASDF 
program package. These files have been consolidated to 
minimize the time required for transfer from DISK02 to 
OSXOOl. The structure of the composite file is listed below: 

BREAKOOl 

EXEC f i le ASDF192 
3REAK002 

TEXT file for POLZRO command 
BREAKOOS 

TEXT file for RESPONSE command 
3REAK004 

Information printed out by the INTRO command 
BREAK005 

VTECl JCL - used for EXEC HRD$CPY 
BREAK006 

VTEC2 JCL - used for EXEC HRD$CPY 
BREAK007 

tables for RESPONSE plots 
BREAK008 

EXEC file to execute the EXEC PUNCH command 
BREAK009 
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Initial file FTOlFOOl with no poles or zeros 
3REAK010 

TEXT file for ,»enerating JOB cards 
8REAK011 

TEXT file for the INTRO command 
BREAK012 

Information printed by OPTIONS command 
8REAK013 

Information printed on the thtle pa^e 
BREAK014 

EXEC file to execute HRDSCPY command 
BREAK015 

TEXT file for the OPTIONS command 
BREAK016 

TEXT file for the softv^are page clear 
BREAK017 
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APPENDIX D. JCL FOR EXEC HRD?:CPY COMMAND 



The Job Control Language required for the execution of 
the EXEC HRD^CPY command is provided to give insight into 
the file manipulation utilized by the ASDF package. Should 
the permanent locations of any of the portions of the ASDF 
package change this section of JCL must be changed also. 



JCL for VTECl 



//GO EXEC PGM=THEPLT,REGION=180K 

//STEPLIB DD UNIT=3330,VOL=SER=niSK03,DISP*SHR,DSN=S1502,ASDF 

//FT05F001 DD DDNAME=SYSIN 

//FT06F001 DD SYSOUT=A 

//FT15F001 DD DDNAME=PLOTPARM 

//PLOTLOG DD SYSOUT=A 

//VECTR2 DD DSN=><i^«VECTR 2, 0 I SP=( , PASS) , SPACE = ( CYL, ( 5, 5) ),UNIT = SYSDA 
//VECTRl DD DSN=^&VECTR1,DI SP=( ,PASS),SPACE = (TRK, ( 1, 1 ) ) , UN I T=SYSDA 
//GO.SYSIN DD * 



JCL for VTEC2 



//PLOT EXEC PGM=I EVMAPP,C0ND=(4,LT) 
//FT15F001 DD DDNAME=PLOTPARM 
//STEPLIB DD DSN=SYS1.VTECPL0T,DI SP=SHR 
//PLOTLOG DD SYS0UT=A 

//VECTRl DD Dl SP=(OLD, DELETE), DSN=^&VECTR1 
//VECTR2 DD Dl SP=(OLD, DELETE) , DSN»&&VECTR2 
//SYSVECTR DD SYS0UT*( A, , 5 55 5 ) 

//VECTTAPE DD DUMMY 

// 
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APPENDIX E: SOURCE DECK FOR ASDF 191 EXEC 
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APPENDIX F: SOURCE DECK FOR ASDF192 EXEC 
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